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THE INITIAL-FINAL MASS RELATIONSHIP: SPECTROSCOPY OF WHITE DWARFS IN NGC 2099 (M37)i 
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ABSTRACT 

We present new observations of very faint white dwarfs (WDs) in the rich open star cluster NGC 
2099 (M37). Following deep, wide field imaging of the cluster using CFHT, we have now obtained 
spectroscopic observations of candidate WDs using both CMOS on Gemini and lris on Keck. Of our 
24 WD candidates (all fainter than V = 22.4), 21 are spectroscopically confirmed to be bona fide WDs, 
4-5 of which arc most likely field objects. Fitting 18 of the 21 WD spectra with model atmospheres, 
we find that most WDs in this cluster are quite massive (0.7-0.9 Mq), as expected given the cluster's 
young age (650 Myr), and hence, high turn-off mass (~ 2.4 Mq). We determine a new initial-final 
mass relationship and almost double the number of existing data points from previous studies. The 
results indicate that stars with initial masses between 2.8 and 3.4 Mq lose 70-75% of their mass 
through stellar evolution. For the first time, we find some evidence of a metallicity dependence on the 
initial-final mass relationship. 

Subject headings: open clusters and associations: individual (NGC 2099) - techniques: spectroscopic 
- white dwarfs 



1. INTRODUCTION 

The initial-final mass relationship connects the mass 
of the final products of stellar evolution for intermedi- 
ate mass stars, white dwarfs (WDs), to their progeni- 
tor mass. It is a required input for the determination of 
the ages and distances of globular clusters from modeling 
their WD cooling sequences (Hansen et al. 2004), for con- 
straining chemical evolution in galaxies, for determining 
supernova rates (van den Bergh & Tammann 1991), and 
for understanding feedback processes and star formation 
in galaxies (e.g., Somerville & Primack 1999). Yet, de- 
spite its fundamental importance, this relation remains 
poorly constrained observationally. 

The first attempt to map the relation was made 
by Weidemann (1977) by comparing theoretical mod- 
els (e.g., Fusi-Pecci & Renzini 1976) of mass loss to the 
masses of a few WDs in both the Pleiades and Hyades 
star clusters. Since then most of the work has focused 
on using observations of WDs in young open star clus- 
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ters to provide empirical constraints on the relationship 
(Koester & Reimers 1981, 1985, 1993, 1996; Reimers 
& Koester 1982, 1989, 1994; Weidemann & Koester 
1983; Weidemann 1987, Jefllries 1997). This effort, span- 
ning almost two decades, is summarized in Weidemann 
(2000). The result is a monotonically increasing relation- 
ship based on about 20 data points, from observations of 
about a dozen star clusters. 

A few recent studies, such as Claver et al. (2001) and 
Williams et al. (2004), have been successful in finding a 
half dozen WDs each in the Praesepe and NGC 2168 clus- 
ters. The synthesis of these data with the earlier studies 
results in an initial-final mass relationship displaying a 
fair amount of scatter. For example, depending on whose 
data one uses, initial mass stars of 3-4 Mq can produce 
WDs ranging anywhere from 0.65-0.8 Mq. This result 
suggests that the relationship may have a stochastic com- 
ponent. 

It is desirable to find a young, rich star cluster with a 
large number of WDs that can be spectroscopically stud- 
ied. NGC 2099 is such a cluster, with 50 WD cluster can- 
didates measured using imaging observations (Kalirai et 
al. 2001a). The cluster is very rich, containing over 4000 
stars, and has a main-sequence turn-off of ^ 2.4 M©. Its 
distance modulus is (rn—M)v = 11.5, and its reddening 
is E{B—V) = 0.23. In this Letter, we present a new WD 
initial-final mass relationship based on 18 WDs in NGC 
2099. In the following Letter, we address the surprising 
result that aU WDs in NGC 2099 (as weU as those found 
in other young open clusters) are all DA spectral type 
(Kalirai et al. 2005). 

2. OBSERVATIONS 

Imaging and spectroscopic observations of NGC 2099 
were obtained with the Canada- France- Hawaii {CFH), 
Gemini North, and Keck I telescopes. In our wide field 
CFHT imaging study (Kalirai et al. 2001a), we found 
~ 67 WD candidates in the central 15' of NGC 2099. 
Based on comparisons between the cluster field and a 
blank field surrounding the cluster, we estimate the rate 
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of field star contamination among our candidates to be 
25% (thus yielding 50 cluster WD candidates). We 
then plotted the locations of all of these faint-blue ob- 
jects on the sky and obtained further observations of 3 
smaller sub-fields in the cluster. These fields were cho- 
sen to maximize the number of WD candidates (all ob- 
jects in the faint-blue end of the color-magnitude diagram 
[CMD] were treated as WD candidates). With Gemini, 
we imaged three 5. '5 x 5. '5 fields using the CMOS multi- 
object imaging/spectroscopic instrument (Murowinski et 
al. 2003). With Keck, we imaged the same three fields 
with the LRIS imaging/spectroscopic instrument that has 
a 5'x 7' field of view (Oke et al. 1995). These imag- 
ing data were not significantly deeper than the original 
CFHT data and were only used to ensure astrometric 
accuracy for the spectroscopy. 

Multi-object spectra were obtained for a single Gemini 
field, and for two of the Keck fields. The Gemini observa- 
tions used the B600 grating, which simultaneously cov- 
ers 2760 (centered at ~ 4700 A). The data were binned 
by a factor of four in the spectral direction to improve 
the signal-to- noise ratio (S/N). The Keck observations 
use the 600/4000 grism (blue side) which simultaneously 
covers 2580 A (centered at 4590 A). For the Gemini field, 
we obtained 22 individual 1-hour exposures spread over 
22 days, taken mostly at low air-masses (< 1.2) and good 
seeing (~ 0."8). For the Keck fields, we obtained 4 2000 s 
exposures in each of the two fields, also at sub-arcsecond 
seeing. 

The Gemini spectroscopic data were reduced using the 
Gemini iraf package, version 1.3. The Keck data were 
reduced using a set of python routines written by D. Kel- 
son (2004, private communication; Kelson 2003). The 
individual exposures were bias subtracted, fiat-fielded, 
cleaned for cosmic rays, wavelength calibrated, sky sub- 
tracted, extracted, combined, and fiux calibrated (using 
bright standard stars) within each of these programs. De- 
tails of the steps involved in each of these procedures will 
be provided in a forthcoming paper (J. S. Kalirai et al. 
2005, in preparation). The only major problem occurred 
for some of the Keck field 2 spectra, which were taken at 
high air-masses and so the bluest flux was lost as a result 
of atmospheric dispersion. Fortunately, two of the three 
stars that we were fitting from this field have also been 
observed in the higher S/N Gemini data. 

In total, we obtained spectroscopy of 24 individual WD 
candidates in the field of NGC 2099 (3 of these stars 
turned out not to be WDs). Therefore, despite sampling 
only 14% of the total cluster area, we include almost 
1 /3 of the total WD population (cluster and field) given 
the careful positioning of the fields. This is therefore 
the largest individual star cluster WD sample that has 
ever been spectroscopically acquired. A CMD showing 
the locations of the NGC 2099 WDs is shown in Fig- 
ure ^ (based on CFHT data). The 18 circles represent 
those objects that we were able to spectroscopically fit 
(see next section). Four of these, shown as open circles, 
represent those objects that have inconsistent theoreti- 
cal magnitudes (from fitting the spectra) as compared 
to the observed magnitudes, assuming they are cluster 
members. Three WDs, which we could not fit with mod- 
els, are shown as crosses. All of the complete spectra are 
shown in the companion Letter (Kalirai et al. 2005). 
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Fig. 1.— Faint-blue end of the CMD of NGC 2099 shown with all 
spectroscopically confirmed WDs in this work (21 objects - based 
on CFHT photometry). The 14 filled circles are those with consis- 
tent theoretical and observed magnitudes (assuming cluster mem- 
bership), the 4 open circles have inconsistent theoretical magni- 
tudes, and the 3 crosses are objects that we could not fit to spec- 
tral models. The model sequences represent WD cooling models 
(Fontaine, Brassard, & Bergeron 2001) at increments of 0.1 Mq, 
with 0.5 Mq at the top and 1.1 Mq at the bottom. 



3. ANALYSIS 

Using the techniques described in Bergeron, Saffer, & 
Liebert (1992), we determine Tcs and log g for each WD. 
The line profiles are first normalized using two points 
on the continuum on either side of each absorption line. 
Therefore, the fit should not be affected by the fiux cal- 
ibration unless there is a strange "kink" or slope change 
at the location of a Balmer line. The fitting of the 
line shapes uses the nonlinear least-squares method of 
Levenberg-Marquardt (Press et al. 1986). The statis- 
tic is calculated and minimized for combinations of Tcff 
and log g, using normalized model line profiles of all ab- 
sorption lines simultaneously. The resulting l-tr errors in 
Toff and log g were tested by simulating synthetic spec- 
tra with the same number of absorption lines, and similar 
S/Ns, and measuring the output parameters from fitting 
these spectra. These results are found to have errors 
slightly less than those in the true spectra (as expected 
given the small errors in fiux calibration and other de- 
fects), and so we use the true spectra errors. Masses 
(Mf) and WD coohng ages (icooi) are found by using 
the updated evolutionary models of Fontaine, Brassard, 
& Bergeron (2001) for thick hydrogen layers ((/(H) = 
Mu/M = 10"") and helium layers of ^(Hc) = 10"^ The 
core is assumed to be a 50/50 C/0 mix. In Figure |21 we 
present the model atmosphere fits for each WD. 

The WD cooling age represents the time that each 
of these stars has spent traversing from the tip of the 
asymptotic giant branch (AGB) down to its present WD 
luminosity. We can now calculate the progenitor main- 
sequence lifetime (tms, the total lifetime of the star up 
to the tip of the AGB) assuming an age for the cluster. 
In Kalirai et al. (2001a), we fitted the NGC 2099 main- 
sequence to solar metallicity isochrones and determined 
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Fig. 2. — Spectroscopic fits (red) shown for the hydrogen Balmer 
fines of 16 WDs in NGC 2099. Within each panel, the lines are 
(bottom), H-y, Hg, H,:, and Hg (top). The number of higher order 
lines varies depending on the spectral coverage of each star. The 
first six objects (WDl - WD7) have the highest S/Ns and were 
obtained using Gemini. The remaining stars (Keck objects) have 
lower S/Ns but have bluer spectral coverage (hence the fitting of 
higher order lines). The identifications are the same as those in 
Kalirai et al. (2005), in which we present the full spectrum for 
all WDs. The corresponding T^ff, log g, and Alf of these WDs 
are given in Table 1. The masses of two additional objects are 
measured as discussed in the text. 



an age of 520 Myr. Recently, C. Deliyannis et al. (2004, 
private communication) have spectroscopically measured 
the cluster metalhcity to be subsolar {Z = 0.011 ± 0.001), 
and the reddening to be E{B-V) = 0.23 ± 0.01, shghtly 
larger than the value that we used. Using these param- 
eters, the age of NGC 2099 is now calculated to be 650 
Myr using the same Ventura et al. (1998) models and 
procedure as described in Kahrai et al. (2001a). We also 
derive similar ages using the Padova group (620 Myr - 
Girardi et al. 2000) and Yale-Yonsei isochrones (630 Myr 
- Yi et al. 2001) for this metalhcity (J. S. Kalirai et al. 
2005, in preparation). 

The tins for each star is determined by subtracting the 
WD coohng ages (tcooi) from the cluster age (650 Myr). 
The initial progenitor masses (Mi) for the WDs are then 
calculated using the Hurley, Pols, & Tout (2000) models 
for Z — 0.01. For this metalhcity, these models give 
very similar values to those derived from the Ventura et 
al. (1998) models. For these masses, a 100 Myr cluster 
age difference would only cause a ~0.3 Mq initial mass 
difference (for 350 Myr lifetimes). This is a very small 
effect on our initial-final mass relationship, and therefore 
the results are not highly sensitive to the derived cluster 
age. Table 1 summarizes the derived parameters for each 
star, with 1-a error bars. Also given are the theoretical 
magnitudes (from fitting the spectra) and the observed 
magnitudes [assuming {m-M)v = 11.5 to NGC 2099]. 
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Fig. 3. — Top - WD initial-final mass relationship shown for the 
18 WDs spectroscopically fitted in this work (circles) and all previ- 
ous constraints (triangles). Also shown are several semi empirical 
and theoretical relations as discussed in the text. Bottom - Closer 
look at those WDs that form the tight sequence with M; = 2.8—3.4 
Mq and are bound by the two Marigo (2001) theoretical relations. 



4. THE INITIAL-FINAL MASS RELATIONSHIP 

In Figure O (top) we present the initial-final mass rela- 
tionship for the sixteen WDs in Figure [21 (filled circles) , 
as well as two other objects, WD6 and WD21 (see Kali- 
rai et al. 2005 for the spectra of these objects). For these 
two stars, we could not determine an accurate log 5 as a 
result of the spectra having too low S/Ns. However, the 
effective temperatures are well constrained, and therefore 
we derive the masses of the stars by combining this infor- 
mation with the luminosities (assuming they are cluster 
members). This gives the radius, which coupled with a 
mass-radius relation, gives the mass of the stars. The 
masses of these two stars, as well as their progenitor 
masses, are found to be in good agreement with others 
in the cluster (WD6 has Mf = 0.92 and Mi = 3.25, and 
WD21 has Mf = 0.85 and Mi = 3.13). Two of the WDs 
shown in Table 1 have negative progenitor lifetimes. This 
is due to the WD cooling age of these stars being larger 
than the cluster age. For WD17, we artificially set its ini- 
tial mass to 7 Mq (the most massive star that produces 
a WD in the Ventura et al. 1998 models). For WD15, 
we compute a 95% confidence lower limit of 408 Myr for 
tcooi and determine the initial mass based on this cooling 
age. These stars are both plotted with open circles and 
an arrow pointing to higher masses to reflect the lower 
limits. It is unlikely that these objects are field WDs 
given their unusually high masses (see Madej, Nalezyty, 
& Althaus 2004 for field WD mass distribution). 

Figure (top) shows that, with just one star cluster, 
we have nearly doubled the number of data points on the 
initial-final mass plane. This is remarkable considering 
the time and effort required to establish the previous 
constraints (shown as triangles). Furthermore, half of 
our data points sit along a very tight sequence with initial 
masses in the range 2.8-3.4 Mq and final masses in the 
range 0.7-0.9 M©. Fitting these data points and their 
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Masses for these two WDs are determined by assuming they are cluster members - see Jl] 
The cooUng time for these two WDs (tcool) is larger than the cluster age - see 



two dimensional errors to a straight line, we determine a 
slope of 0.33 and an intercept of -0.19. The residual rms 
of the fit is 0.03. This suggests that stars with masses 
between 2.8 and 3.4 Mq, and metallicities slightly less 
than solar, will lose 70-75% of their mass through stellar 
evolution. 

Four of our stars are seen well below the tight sequence 
of points discussed above. Although only one of these 
four stars has an inconsistent theoretical magnitude, they 
could still be field WDs located near the cluster. Not 
only are their masses consistent with the distribution of 
field WDs (M = 0.56 Mq - Sloan Digital Sky Survey 
data; Madej, Nalezyty, & Althaus 2004), but we also 
expect 4-5 WDs in our sample to be field stars based 
on WD number counts in blank fields surrounding NGC 
2099 (Kalirai et al. 2001a). The lower masses of these 
stars could also be attributed to binary star evolution 
and mass transfer. We will investigate these possibilities 
further in J. S. Kalirai et al. (2005, in preparation). 

We have plotted several theoretical and semi empirical 
relations in Figure |3| The dotted line is the revised semi 
empirical initial-final mass relationship from Weidemann 

(2000) . The dashed line is the theoretical initial-final 
mass relationship from the Padova group stellar evolu- 
tion models (Girardi et al. 2000; L. Girardi 2004, pri- 
vate communication) for Z = 0.008. The two solid lines 
are theoretical initial-final mass relations from Marigo 

(2001) , for Z = 0.008 (top) and Z = 0.02 (bottom). The 
latter models have improved mass-loss mechanisms in 
their post-main sequence evolutionary phases and there- 
fore should be preferred over the others (L. Girardi 2004, 
private communication). Considering this, it is inter- 
esting to note that half of our NGC 2099 cluster WDs 
reside in a region of the initial-final mass relationship 
bound by the two Marigo (2001) relations (see bottom 



panel of Figure |3J|. This is what we would expect, given 
that the cluster metallicity falls between that used in the 
two models {Z = 0.011). We point out, however, that 
the errors in our measurements are too large for us to be 
able to discriminate between the two relations. Our data 
points are, however, mostly located above the other pre- 
vious constraints. This can also be explained by the fact 
that almost all of the other clusters are solar or higher 
metallicity, resulting in more efficient mass loss during 
stellar evolution (Marigo 2001). Although this could be 
a systematic effect (e.g., the masses of the WDs in other 
clusters have been derived using different models), these 
results may be suggesting that, for the first time, we are 
seeing the effects of metallicity on the initial-final mass 
relationship. In J. S. Kalirai et al. (2005, in preparation), 
we will eliminate these systematics by reevaluating all of 
the previous data on the initial-final mass relation (in- 
cluding the new results of Williams, Bolte, & Koester 
2004) using a consistent fitting method and set of evolu- 
tionary models. 

The spectroscopic measurement of masses of WDs with 
y ~ 23 has not been previously accomplished. The num- 
ber of targets that are accessible at these magnitudes is 
several orders of magnitudes larger than previously iden- 
tified (see work on the on-going CFHT Open Star Cluster 
Survey; Kalirai et al. 2001a; 2001b; 2001c; 2003). With 
continuing observations of rich open star clusters, we can 
envision placing >100 data points into this very funda- 
mental relation. By observing both younger and older 
clusters, the entire initial progenitor mass range can be 
constrained and a detailed initial-final mass relation can 
be produced. 
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